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1 Introduction 

The string/gravity correspondence allows to calculate correlators in certain strongly coupled gauge 
theories via solving the equations of motion of supergravity (SUGRA). Similar to calculations in 
field theory, this requires regularization and renormalization. In the case of asymptotically AdS- 
spaces, the method of holographic renormalization (HR) has been developed systematically in 
[Il[2l[3l[ll[5l[6l[71[8l[9]. However, these methods do not cover cases in which the field theory has a 
logarithmically running coupling even in the ultraviolet (UV)§ Holographically, this translates to 
a bulk metric which is not asymptotically AdS (aAdS), but has logarithmic warping (in a suitably 
defined radial variable). The prime example of such a background is the Klebanov-Strassler (KS) 
solution [11], which is well approximated in the UV by the Klebanov-Tseytlin (KT) solution |12j . 
Calculating correlation functions in these cases is much more involved, also because the procedure 
of HR has not been worked out yet in a systematic way similar to the aAdS case. 

^Recently, progress was made on the holographic renormalization in bulk backgrounds which are conformal to 
AdSp+2 X S^~'' with a non- vanishing dilaton [10] , These cases imply couplings that run with a power law in the UV. 
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As a result, only a few attempts to calculate correlators using the KT background have been 
made, cf. [131 El [ISl [l6], and only in [15l HB] the program of HR, as reviewed in [7], was applied. 
Furthermore, calculations of mass spectra in the KS background |131 \T7\ [T8| [T9l [20l [21] have been 
done using a pragmatic approach assuming that a consistent method of HR in aAdS backgrounds 
exists. 

In this note, we would like to readdress the question of how to calculate renormalized correlators 
holographically from backgrounds which are not aAdS, given that this is a feature one would expect 
for the dual description of any gauge theory with a running coupling in the UV, like QCD. Our 
approach is different in spirit from [151 [16] . We consider a general bulk theory of gravity coupled 
to an arbitrary number of scalars, whose potential can be expressed via a "superpotential" . Such 
theories are known as "fake SUGRA" theories [22^ and allow for BPS domain wall background 
solutions, which are the holographic duals of renormalization group flows. The fake SUGRA systems 
include the case of KS (and also KT), when viewed as a consistent truncation of type-IIB SUGRA 
|25[ 126] ■ For such a general theory, it would be a daunting task to find the complete counterterms. 
Thus, we take a step back and content ourselves with giving a recipe how to calculate renormalized 
two- (and to some extent one-) point functions of the operators dual to the scalars of the theory. 
Furthermore, we only consider field theories living on a flat space-time, which allows us to ignore 
all counterterms involving the space-time curvature. In a sense, our approach is inspired by [8l[9], 
where the philosophy was put forward to concentrate on the part of the counterterm action which 
is really necessary to calculate ?i-point functions for a given n, i.e., the terms of n-th order in the 
fluctuations. In this spirit, we consider the case n = 2. The counterterms we propose involve the 
fluctuations in a covariant way, but, otherwise, do depend on the background. It might be possible 
to derive them from a fully covariant expression, but we have not attempted to do so. 

The starting point of the holographic calculation of correlation functions in AdS/CFT is the 
correspondence formula [27] 



where S'on-shls] denotes the renormalized bulk on-shell action evaluated as a functional of suitably 
defined boundary values Sj of the various bulk fields, which are identified with the sources coupling 
to certain QFT operators Oi (we will make this more precise below). Hence, the bulk quantity 
'S'on-sh[s] is identified with the generating functional of the connected correlation functions of various 
QFT operators. In particular, the exact one-point functions of the QFT operators are given by 







(1.2) 



^Here "fake" does not mean that the theory is necessarily non-supersymmetric, just that the formalism is applicable 
more generally. The relation between supergravity and fake supergravity was analyzed in [231 124j . 



2 



In order to calculate the two-point functions, one has to know the dependence of the right hand 
side of ()1.2p on the sources sj up to linear order. Determining this dependence requires to solve 
the linearized bulk equations of motion. Thus, if we are interested in one- and two-point functions, 
a knowledge of the action, which consists of bulk and boundary terms, up to quadratic order in 
fluctuations is sufficient. We will give a recipe how to calculate the quadratic terms in section [3?n In 
doing so, we make use of the gauge-invariant formalism for the fluctuations developed in [28^ 129 1 [26]. 
in which the scalar fluctuations explicitly decouple from those of the metric at the linearized level. 
Thereby, the gauge-invariant fields are identified with the relevant bulk degrees of freedom that 
encode the information on the boundary correlation functions. We restrict our attention to the 
scalar sector, but the recipe can be extended easily to the traceless transversal fluctuations of the 
metric. 

In order to calculate the vacuum expactation values (VEVs), one would also have to know the 
boundary terms linear in the fluctuations |f| We do not know yet how to generalize our prescription 
to those terms. Thus, strictly speaking, we can only calculate the contributions to the one-point 
functions, which are linear in fluctuations (i.e. excluding the VEVs). However, we shall observe 
that the linearized equations of motion have a zero mode solution (depending only on the radial 
coordinate) which seems to encode some information about the VEVs. We will make this more 
explicit in the examples that we discuss in later sections. 

One advantage of our approach is that it allows to discuss the scheme dependence of one- and 
two-point functions in a rather general way, as we will do in section 13.31 and then more concretely 
in the examples. 

The outline of the rest of the paper is as follows: In section [2l we review the linearized bulk 
dynamics of a system of scalars coupled to gravity, expressed in the gauge-invariant variabels of 
[28\ [29l 126] . Section [3] contains our main results. It gives the general recipe of how to calculate 
renormalized two-point functions and makes some comments on the VEV part of the one-point 
functions. It also discusses the issue of scheme dependence. As a first check, we then continue 
in section d] by applying the general formulas of section [3] to two aAdS backgrounds, which have 
been studied extensively in the literature. These are the GPPZ [3D] and the Coulomb branch flows 
[3H 132]. In these cases, we find complete agreement with earlier results, including the full scheme 
dependence of the two-point functions. This gives us enough confidence to carry on and consider 
the KS case in section [5j The analysis requires lengthy expressions for the asymptotic solutions of 
the coupled linearized scalar equations of motion, which we include in appendix [Al 

■^Of course, also a term independent of the fluctuations has to be added in order to obtain a finite action. Through- 
out the paper we assume that such a term has been added. 
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2 Bulk Dynamics 



Let us start by reviewing the equations governing the dynamics of the bulk fields [291 [26] , which 
encode the information about two-point functions in holographic renormalization group flows. 
The systems we consider are fake SUGRAs in d + 1 dimensions with actions of the form 



S 



id+l 



-^R + ^Gab g 



MN 



(2.1) 



with M, N = 0, 1, ... d, and where the potential V{(l)) is given in terms of a superpotential W^cp) 
by 



1 



d 



(2.2) 



We will not specify at this point the boundary terms in ()2.ip . as they do not affect the bulk 
dynamics, although they are important for holographic renormalization. Our notation agrees mostly 
with [2^. In particular, field indices are covariantly lowered and raised with the sigma-model metric 
Gab and its inverse, G"*, respectively, Wa = daW = dW {(/))/ d (p"" , and covariant derivatives with 
respect to the fields are indicated by Da or by a "|" preceding the index, as in Wa\b = D^Wa = 
dbWa — G^ab^c, Q'^ba being the Christoffel symbol for the metric Gab- 

Holographic renormalization group flows are described by domain wall backgrounds of the form 



ds2 = dr2 + e2^(^) rif,^ dx^ dx^ , 4>'' = ^'^(r) , 



(2.3) 



with fijU = 1, . . . d, satisfying the BPS equations 

2 



drA 



d-l 



(2.4) 



Linearized fluctuations around such a domain wall background are best described in a gauge in- 
variant fashion, in which thi 
and the scalar fluctuations 



variant fashion, in which the independent fields are the traceless transversal metric fluctuations, e*-, 



h 



(2.5) 



Here, are scalar field fluctuations and h is the trace part of the metric fluctuations. Furthermore, 
0{f'^) denotes corrections which are non- linear in the fluctuations. The fields a" and e*- satisfy the 
(linearized) equations of motion 



and 



Dr + M 



dr 



2d 



W ] [Dr - M) + e-2^ □ 



= 



2d 



W]dr + e 



~2A 



□ 



e^- = 



(2.6) 



(2.7) 
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respectively. In (j2.6p . we have omitted the field indices, M denotes the matrix 



W 



(2.8) 



W 



and Dr is the background covariant derivative 



DrO! 



a 



(2.9) 



For more details, we refer the reader to the original papers [28 \ 129^ [26]. 

In this paper, we focus on the scalar field equation (j2.6p . Let rig be the number of scalar fields 
(components of a). As in [19], we shall assume the existence of a set of 2ns independent solutions of 
(j2.6p . which are defined as power series in k'^ (in momentum space), with r-dependent coefficients 
that are more and more suppressed with increasing powers of fc^Q Moreover, the leading term (for 
large r) in each solution should be independent of k'^. In position space, k'^ simply translates to the 
operator — □. Amongst these solutions, one can distinguish between Ug asymptotically dominant 
solutions O-i {i = 1, . . . ,ns) and ng sub-dominant solutions with respect to their behaviour at 
large r. Including the field index, we shall interpret df and d? as Ug x Ug matrices. A regularity 
condition in the bulk interior allows only for Ug independent regular combinations of the asymptotic 
basis solutions. Hence, we shall decompose a general regular solution of (|2.6p into 



where Sj and tj are called the source and response coefficients, respectively, and = V^'^'^'^- 
The bulk regularity condition uniquely determines the (functional) dependence of the responses tj 
on the sources Sj and gives rise to the non-local information for the two-point functions of the dual 
operators. 

Throughout the paper, we shall consider mostly the analogue of (j2.10p in momentum space, 
sometimes omitting the dependence on k. Moreover, a • will be used to denote the inner product 
in field space, or the contraction of field space indices, e.g., a • b = a°'Gabb^- 

3 Perturbative Holographic Renormalization 
3.1 Scalar Two-Point Functions 

In this section, we shall present the general formalism for obtaining finite, renormalized two-point 
functions for the QFT operators that are dual to the bulk scalar fields. Our starting point is the 

^This is tantamount to demanding that the warp function A(r) grows without limit for r oo, so that e^^^ 
in (|2.6p can be regarded as a correction in the asymptotic region. 



a"(r, x) = d^(r, -□^)5i(x) + d^(r, -D^) Xi{x) 



(2.10) 
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following action, which is quadratic in the fluctuations and encodes the bulk field equations (j2.6p . 

5 = i y d'^+^x e'^^iUDr - M)a] ■ [{Dr - M)a] + e'^^^^a • a^o} + d^x e*^^ a • C/ • a , (3.1) 

with some symmetric counterterm matrix U, which is a local operator that will be specified in a 
moment. The bulk integral in p.ip is to be understood with a cut-off tq, where also the boundary 
counterterm is evaluated. It follows that the variation of the on-shell action with respect to a 
variation of the boundary value o"(ro) is given by 

^-^ = e'^{Dr-M + U)aa, (3.2) 

where the right-hand side is evaluated at r = tq. Let us define the counterterm matrix as 

Uab = Mab-^ [{DraU{a-%b + {Dra)ib{a-^)ia] , (3.3) 

where {a~^)ia is the inverse of the matrix af, defined in momentum space as a series in k'^, or 
equivalently, in position space as a series in — □. We will see momentarily that this definition 
leads to finite one- and two-point functions. We also note the following subtlety. The counterterm 
in p.ip needs to be local in the fields, which means that Uab should be a polynomial in k'^ (in 
momentum space) or — □ (in position space). The assumptions made in section [2] imply that Uab is 
a series in k'^. However, we also assumed that the coefficients of the series a with increasing powers 
of k"^ are suppressed for large r due to the factor e"^"^^'"^ so that we can truncate the series in (|3.3p 
to some polynomial, because the terms thus neglected vanish in the large-r limit. Hence, strictly 
speaking, the counterterm operator Uab in (|3.ip is a polynomial truncation of ()3.3p . 
Before deriving the two-point function, let us also introduce the following matrices 



(3.4) 



These matrices are independent of r, as one can show from the field equation (j2.6p . This implies that 
Zij should be identically zero, as the sub-dominant solutions vanish fast enough asymptoticallyl^ 
Furthermore, they are functions of k"^ or — □, depending on whether one works in momentum or 
position space. 

Combining (j3.2p with the decomposition (j2.10p . the (linear term of) the exact one-point function 



Z.J = e''^ 


[{Dra)i ■ 




- ai ' 




Z,j = e'^ 


[iDra)i ■ 




- ai ' 


■ {Dra)j 




[{Dra)i ■ 




- ai ' 





^This is not necessarily the case if there are two or more bulk scalars with mass = 2(2 — d), which, in the 
aAdS-case, would be dual to operators of dimension A = 2. If at least two of these scalar fields are present and the 
background is not aAdS, one has to check more carefully whether Zij indeed vanishes. We will assume this in the 
following, as it simplifies our final result. In all the examples we are considering later, this issue does not play any 
role. 
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1.2p . in momentum space, takes the for 



nfl 



(^^^(A;))^ = - lim e'^^^'^ 



{Dr- M + U)[ai5i{k) + aiXi{k)] , (3.5) 



where, for the sake of brevity, we have omitted the dependence of the asymptotic solutions o and 
d on r and k"^. 

Substituting ()3.3p into (j3.5|) and using the matrices (j3.4p . after some algebra one obtains 



1 ~ 1 dx 

{Oi{k))^ = Z.jXj + -Zij5j + ^Zjkg^Xk 

1 / dxk dx- dxk \ ^"^'^^ 

Here, we have omitted the arguments k on the right hand side. To obtain the final result, we observe 
that the third term on the right hand side vanishes, since Zij = 0, as stated above. Moreover, the 
last term, which is the only one with a cut-off dependence, vanishes when the large-r limit is taken, 
because (a~^); • dfc goes to zero. Hence, we end up with 

(C'i)i = ZijXj + -ZijSj , (3.7) 



which holds both in momentum and position space. From (j3.7p . one obtains the connected two- 
point function 

{Oi{x)0,{y)) = Za.(-n.)|^ + ^4(-n.)<5(x - y) . (3.8) 

As promised, (13. 7p and (13. 8p are finite in the limit rg — > oo, as the matrices Zij and Zij do not 
depend on r. Eq. (j3.7p agrees with eq. (2.24) of [19], for which it provides the missing piece Yij = Zij 
and identifies the contact term, which shall be discussed further in subsection 13.31 

In momentum space, ()3.8p has a more practical form. Setting y = by translational invariance 
and Fourier transforming the coordinate x, one finds 

I d'^xe^^^ {0,ix)0,m = Z,k{k^)^{k) + ^Z,,{e) . (3.9) 

In what follows, we shall often work in momentum space omitting the argument k. By the two-point 
function (OiOj) in momentum space, we will intend (j3.9p . 

Unfortunately, the symmetry of the two-point function under exchange of Oi and Oj is not 
obvious from (j3.8p . It would be a non-trivial test of any concrete calculation to see whether the 
right hand side, with its antisymmetric second term, combines to something symmetric. 



®The subscript 1 on the left hand side indicates that these are just the terms linear in the fluctuations. 
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3.2 VEVs 

Let us make a few comments on VEVs. Equation (j3.5p only gives the part of the one-point function 
which is hnear in the fluctuations. At the moment, our approach does not allow for a systematic 
derivation of the VEVs yet. However, we would like to make the following observation. The scalar 
equations (|2.6p have the zero mode solution 

a"^ = G-^"^ , (3.10) 

which only depends on the radial variable r. Like any fluctuation, this has an expansion as in ()2.10p 
(with vanishing k"^). We will see in section H] that the response function r of this zero mode solution 
encodes the VEV for the GPPZ and the CB flows. It would be interesting to understand how 
generally this holds. In section [5] we will discuss the issue of VEVs in the KS theory by analyzing 
the response functions of the corresponding zero mode solution. 

3.3 Scheme Dependence 

In QFT, contact terms of correlation functions, which do not influence physical scattering ampli- 
tudes, depend on the renormalization scheme. Let us now discuss how the scheme dependence of 
the two-point functions (jS.Sp appears from a bulk point of view. For the sake of brevity, we shall 
work in momentum space and omit all functional arguments. 

The starting point is the decomposition (|2.1Up of a regular solution to the bulk field equations. 
Clearly, the definition of the asymptotic solutions dj and dj is not unique. Our restriction on 
the functional form of these solutions in terms of series of k'^ and the fact that all sub-dominant 
solutions are negligible for large r with respect to all dominant ones still allows for a change of 
basis of the form 

d- = AijO-j + XijO-j , d- = iiijUj , (3-11) 

where the (non-degenerate) matrices Aj^, Ajj and fiij are polynomials in /c^. Under this change of 
basis, the matrices Zij and Zij transform intcQ 



^ij — Kk^jlZkl + {Kk^jl - ^jkKl) Zkl , 

respectively, while the source and response coefficients in the new basis become 



(3.12) 



s'i = 5j{A-^)ji , x'i = [tj - 5i{A-^)ik\kj] ■ (3.13) 

Inserting these transformations into (|3.8p . one finds the connected two-point functions of the oper- 



^Remember Zij = 0. 
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ators 0[ coupling to the sources s[, 



{p',0'>) = K,k^ji {OkOi) - - {AikXji + AjkXii) Zki . 



(3.14) 



Hence, the matrix Aij performs a rotation of the basis of operators, as one would have expected, 
while a non-zero Xij changes the contact terms, which corresponds to a change of renormalization 
scheme. 

In QFT, operators are usually characterized by their scaling dimension, which is renormalization 
scale dependent. Under renormalization, they undergo operator mixing, such that an operator of a 
given dimension, defined at a certain renormalization scale, is generically made up of the operators 
of equal and lower dimensions, defined at a larger renormalization scale. There is, however, some 
ambiguity, as operators of equal dimension and otherwise equal quantum numbers can be arbitrarily 
combined to equivalent combinations. This ambiguity finds a natural counterpart in the present 
approach. Ordering the dominant asymptotic solutions according to their asymptotic behaviour in 
descending order, it is natural to choose A in "upper triangular" form, such that each dominant 
solution is modified only by solutions of equal and weaker asymptotic behaviour. A similar remark 
would apply for /i. 

A further restriction on the redefinition could come from the fact that the lowest order terms in 
a near boundary expansion of the dominant solutions typically have a definite correlation between 
powers of e~'^ and powers of k^. This is well known, for instance, in the aAdS case with a single 
scalar field, cf. the discussion in sec. 5.1 of [7]. Something similar is apparent from the asymptotic 
solutions in the case of KS, cf. appendix lA.li We shall refer to a choice of dominant solutions 
respecting this correlation "natural" choice. 

Finally, we remark that it is reasonable to assume that Aij and/or /ijj can be chosen such that 
Z'-j = 6ij in (j3.12p . A possible obstruction to this possibility would be that the matrices needed to 
achieve that are non-polynomial in k'^. We shall see later that the choice Zij = 6ij is possible for 
the KS system. Starting with such a choice, a further change of basis using just Xij would lead to 



implying that one can achieve Z'^j = by a suitable choice of Xij , although this choice is obviously 
not unique. 

4 AAdS Examples 

In this section, we shall compare our results from the last section with the results of holographic 
renormalization in the case of asymptotically AdS bulk space-times. Our reference for holographic 
renormalization will be [6], as it is the most complete concerning scheme dependence. In particular, 




(3.15) 
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we shall consider the active scalars in the GPPZ [30] and the Coulomb branch flows \31\ I32j. 
4.1 GPPZ Flow 

The GPPZ flow [30] (d = 4) is characterized by the superpotential 

T^(<^) = -^(cosh^ + l) . (4.1) 
The bulk scalar fleld (j) is dual to an operator of scaling dimension 3 and behaves asymptotically as 

cj){r) = (t>o e-" +^2r e"^" +(/)2 e'^^ + • • • , (4.2) 
where (t)Q and 4)2 are independent coeflicients, and -02 = ^n</'o- The background solution satisfles 

(4.3) 



1 -e 
which implies 



T 

aim- 



00 = ^/3 , 02 = , ^2 = ^- (4.4) 
le exact one-point function of the operator O coupling to the source 0o is given by eq. (72) 

{O) = 202 - (mo + ^ 000 + y 0i] , (4.5) 



where the two scheme-dependent coefficients mo and M4 stem from the addition of the finite coun- 
terterms 

j d4xV^(^04 + imo5^'^a^09.0) . (4.6) 
Renormalization schemes with = —4:/2> respect SUSY, but ttiq remains undetermined by SUSY 

Let us now consider an arbitrary scheme. Linearizing around the background (|4.3|) . = + y?, 
and switching to momentum space, the exact one-point function (j4.5p becomes 



2 V 3 



3 / 4\ / 1\ 



2(932 - V^o) + 2 ( ^4 + 2 1 V^o + ( "^0 + 2 ) ^ V^o 



(4.7) 



The first term on the right-hand side is a scheme-dependent VEV, which vanishes in SUSY schemes. 
The term in the square bracket encodes the 2-point function. 

To proceed, we need to relate the field fiuctuation to the gauge-invariant variable a. The 



*We have omitted the curvature-dependent terms, which are irrelevant here. 
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necessary details can be found in [26]. In particular, eq. (4.19) of that paper and (j2.5p read 

where W = dW/ d(j) and z/ = in the orthonormal gauge, which was used to derive (j4.5p . Moreover, 
the fields b and a are related on-shell by (4.41) of [26| . 



w 

b = --a . (4.9) 
Using these relations, one can derive the leading behaviour of the gauge-invariant field a as 

a(r) = (^oe-''-^A:Vor-e-3^+(992 - ao)e-3^ + --- , (4.10) 
from which one reads off tpo = ao and (/?2 = 12 + Oo- Hence, ()4.7p becomes 

2,-^.3,. 202 + ^ + ^ j ao + f ?no + ^ j /c^oo . (4.11) 



Let us now compare this expression with the results of section [3711 In particular, let us define 
the dominant and sub-dominant solutions by 

a = e-^-re-3^^A:2 + e-3^a2 + -- - , d = e'^^ + • • • , (4.12) 

with an as yet undetermined coefficient Q2- From these, one obtains the counterterm "matrix" 

U = -r e-2^- k'^ + e-2'^ (-k'^ + 202] +■■■ , (4.13) 



2 

where the ellipses indicate terms which do not contribute and can be truncated. It can be noted 
that the first term on the right-hand side agrees with the standard logarithmic counterterm [6] . The 
terms in brackets give finite contributions. Moreover, ()3.4p yields Z = and Z = 2. Expressing 
(j4.10p in the basis (j4.12p . the source and response coefficients become 

s = ao , t = a2 - a2ao , (4.14) 

respectively, so that (|3.7p yields 

(O)i = 2(a2-a2ao) . (4.15) 

Comparing this with the linear term in (|4.1ip . we find agreement in the non-local term containing 
a2 and determine 02 as 

3 / 4\ 1,„ ..,2 



a2 = -^U4 + -j --(2mo + l)/c' . (4.16) 
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This result states explicitly the relation between the choice of the dominant basis and the QFT 
renormalization scheme. Note that only 02 ~ (/c^)" for n = or 1 lead to a change of scheme. 
Extending our discussion in section 13.31 we might call these /c^-dependences for redefining the 
dominant solution also "natural" choices. In particular, a2 = in the SUSY scheme {u^ = —4/3) 
with mo = —1/2, which agrees with (5.23) of [7]. 

Finally, let us check that the VEV is encoded in the response coefficient of the zero-mode 
a = W'/W = ^e-\ This imphes % = 2/^3 and a2 = 0, and from (fiTi]) and (f4T6]) follows (here, 
F = 0) 

2 ^/3 f 4' 



r = -^a2 = ^^U4 + -j . (4.17) 

The right-hand side coincides with the constant term in (j4.1ip . i.e., in any renormalization scheme, 
the response coefficient of the background mode is just the VEV. In the SUSY schemes, the VEV 
vanishes, as is well known for the GPPZ flow. 

4.2 Coulomb Branch Flow 

Another interesting example is the supergravity dual of a particular state in the Coulomb branch 
of the = 4 SYM theory [311 [32]. In this case, the superpotential is given by 

_ 20 1 ■4'#> 3 „ „ 

VF((/,) = -e ^ --eV6 = __ -</,2 + 0(03) _ (4^8) 

The corresponding potential can be easily computed using ()2.2p resulting in 

y(<^) = _3 _ + J_^s ^ 0(04) _ (4_19) 

Thus, we see that (/> is dual to an operator of scaling dimension A = 2, as it has a mass = —4. 
Correspondingly, it has an asymptotic expansion 

</.(r) = (/)ore-2^+(^oe-2- + ... , (4.20) 

where (po and <j)o are the two independent coefficients. The background solution is given by 

We remind the reader that this domain wall solution can be lifted to 10 dimensions, where it 
describes a uniform distribution of D3-branes over a disc of radius /. At the boundary, the scalar 
field vanishes at the rate 

0«--^^%-^^ (4.22) 
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implying 

P 

<Po = 0, (1)0 = -^. (4.23) 
The exact one-point function of the corresponding operator is given by [33^ [U [6] 

{O) = 00 + ^2(^0 • (4.24) 

Two comments are in order here. The second term, involving a scheme dependent constant U2, is 
new compared to the corresponding formulas of [331 El E] and arises from the addition of a finite 
counterterm proportional to 

$ . (4.25) 

r=ro 

Furthermore, our result differs from [33\ [5] by a factor of 1/2, as our (j)Q differs by a factor —2 from 
theirs and our definition of the one-point function (|1.2p exhibits an additional minus sign. 
Linearizing around the background, (p = (f) + ip leads to 

(O) = -— + ^o + U2^o . (4.26) 

The first term on the right hand side constitutes a finite VEV of the dual A = 2 operator, which 
is independent of the renormalization scheme. In order to express the first-order terms in (j4.26p 
in terms of the gauge invariant variables, we make use of the relations (j4.8p and (j4.9p again. This 
results in 

a = ip + 0{e-^'') . (4.27) 

Using the definition 

a = aore-^^+aoe-^^ + ... , (4.28) 
the one-point function in gauge invariant variables reads 

P 

{0) = --^ + ~ao + U2ao . (4.29) 
In order to make contact to section [3Tl let us consider the dominant and sub-dominant solutions 

a = r e-^' +a e''^' , a = e'^^ . (4.30) 
Writing (I4.28P in this basis, one can read off the source and response as 

s = ao , r = ao — aao . (4-31) 
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Moreover, in this case (j3.4p results in Z = and Z = \. Thus, (j3.7p reads 

(O)i = ao-aao. (4.32) 

Comparison with the part of (|4.29|) which is hnear in fluctuations implies that a = —u^. 
According to (j3.3p the counterterm "matrix" obtained from the basis (j4.30p is 

C/ = -- + 4 + C^(^"') • (4-33) 

This provides the standard logarithmically divergent counterterm [6] and a scheme dependent finite 
contribution. 

Finally, let us again consider the zero-mode a = In this case, it is given by 

W' 4 /2 

Hence, we get ao = ^i^d Oq = 0, and thus, the response is 

4 P 

Comparing with ()4.26p . we see that, in the CB flow, the response of the zero-mode gives the VEV 
only up to an overall factor, but it is non-vanishing independently of the scheme. 

5 KS System 

Finally, we would like to apply the general discussion of section [3] to the case of the Klebanov- 
Strassler theory. We start by reviewing the relevant facts about the background, but we will be 
brief and refer for more details to |26^ [T9] and references therein. 



5.1 KS Background 

The effective 5-d model describing the bulk dynamics of the KS system contains seven scalar fields 
We shall use the Papadopoulos-Tseytlin [25] variables ?/, 6, /ii, /i2)- The dual operators have 
dimensions A = 8, 7, 6 and twice 4 and 3 each0 







®As in [261119) . we restrict ourselves to the J^'~' = 0^^ scalar sector, where C denotes the quantum number under 
the Z2 charge conjugation symmetry of the KS theory, cf. [341 135) . Additional scalar fluctuations with J^'^ = 0+" 
and = — were discussed in [54ll35][20] . 

^"it is not completely clear to us whether this is a meaningful statement, because in contrast to aAdS settings, the 
KS system has no UV conformal flxed point, where the operator dimensions can be fixed. However, the deviation 
from aAdS behaviour is quite mild, such that the asymptotic solutions behave nearly as if the dual operators had 
definite dimensions. This can be seen explicitly by inspecting the asymptotic solutions given in appendix [X] Their 
exponential r-dependence (e''^"*''^''^ for a and e"^'^''^ for a) is what one would expect for a solution dual to an 



14 



The sigma-model metric is 



+ \ e-*-'" [e-=^^ dM{h, - h2)d''{h - h2) + dM{hi + h2)d^\hi + /12)] , (5.1) 

and the superpotential reads 

^^ = -^(e-2f-2- + e4fcosh2/) + ie^P-2"[g + 2P(6/i2 + /ii)] . (5.2) 

Here, Q and P are constants related to the number of D3-branes and wrapped D5-branes, respec- 
tively. 

It is useful to introduce the KS radial variable r by 

dr = e^P dr . (5.3) 

In terms of r, the KS background solution of ()2.4p is given by 

^> = $0 , (5.4) 
ey = tanh(T/2) , (5.5) 

sinhr ' ^ ^ 

/ii = + Pe*Ocothr(Tcothr- 1) , (5.7) 

„ <f> T coth T — 1 , , 

/i2 = Pe*« — , (5.8) 

smhr 

2 g6p+2x ^ ^^^^ ^ ^ ^ ^5 

3 sinh r 



with 



00 



h^r) = [ M ^ [2 sinh(2^) - 4^?]^/^ ^ ^^^^^^ 

y sinh 

T 

Moreover, the warp factor is given by 

^ e^*' (e-2- sinhr)'/^ /i(t) , (5.12) 
with a proportionality factor that sets the momentum scale. 

operator of dimension A [the KS radial variable r will be introduced momentarily in (|5.3|l ]. Thus, we still regard the 
concept of dimension as useful for distinguishing the different asymptotic solutions and we will employ the expression 
correspondingly. 
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The KT background solution is somewhat simpler, because there ?/ = 6 = /12 = 0, but it has a 
singularity. For the KT background solutions of the other fields, we refer to |26j . 

The gauge- invariant fluctuations around the KS background fulfill the field equations (j2.6p . 
whose explicit form can be found in [19]. As in that reference, we form the sigma model covariant 
fluctuations a'* using ()2.5|) for the rescaled fields 



In this way, there is no explicit dependence of the matrix M [defined in (j2.8p ] on P and The 
resulting equations couple all the seven scalars non-trivially, but can be solved iteratively for large 
r. The asymptotic solutions are collected in appendix lAl 

The astute reader will notice a pattern in the solutions of the appendix. Let us consider the 
two groups of scalars consisting, on the one hand, of hi and and on the other hand, of y,b 
and /12CJ In |19j . these two sets of scalars were called the "glueball sector" and the "gluinoball 
sector", respectively. In the KT background, the scalars in the gluinoball sector are inert, i.e., their 
background solutions are identically zero, and consequently any terms coupling the two sectors are 
absent. This eventually leads to the singularity in the IR, which is resolved in the KS background 
by taking into account the backreaction on the gluinoball sector. Nevertheless, the UV decoupling 
is also apparent in the asymptotic solutions of the appendix. The dominant solutions 01,03, 04 and 
05 and the subdominant solutions di, 03, 014 and 05, which are related to the operators of dimensions 
A = 8, 6 and 4, only have the first four components excited at leading (and next-to-leading) order. 
These four components correspond exactly to the scalars of the glueball sector. The mixing only 
appears at order e^"^ relative to the leading order, as is to be expected from the asymptotic expansion 
of the equations of motion, cf. section 5.4 in [TUj. Similarly, the dominant solutions 02,06 and 07 
and the subdominant solutions 02, dg and 07, which are related to operators of dimensions A = 7 
and 3, only have the last three components excited at leading (and next-to-leading) order. These 
correspond to the scalars in the gluinoball sector. 

5.2 Holographic Renormalization 

We are now ready to apply the formalism of section [3] to the case of the KS system. In the following 
discussion, we often restrict ourselves to the A < 4 operators. This simplifies the calculations 
considerably, and it suffices to discuss all the general features of a system with several coupled 
scalars. 

In order to discuss the issue of scheme dependence, we allow for redefinitions of the dominant 
asymptotic solutions given in appendix lA.il with the subdominant solutions of appendix I A. 2 1 Again, 
^^More precisely, we consider the gauge invariant scalars built on them according to (I2.5|l . 




(5.13) 
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we restrict to the A < 4 operators and only modify the corresponding dominant solutions according 
to 



dj + Ajod 



(5.14) 



with i,j = 4, 5, 6, 7. 

The matrix Z from ()3.4p can be calculated using the solutions of appendix |A] and is given by 
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4 
9 

















2 
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5 
9 














1 ) 



(5.15) 



where we have introduced the abbreviation 



3V3 

m 



(5.16) 



Note that Z does not depend on the Xij, i.e., it is scheme independent as it should be according to 

(nmi) . 

As discussed in section [3]3l one can also redefine the dominant solutions by other dominant ones. 
In particular, to a dominant solution of dimension A, one could add other dominant solutions of 
dimensions smaller than or equal to A. This would amount to an upper triangular matrix A, cf. 
(j3.1ip . It is easy to check that using fiij = 6ij and 



A,,- = f3V3p4g2</,o 



-1 



/-A 

80 










V 





_ 9 
2 









1280^ 


59697 p2 
204800 ^ 


3051207 fj2 
2048000 ^ 





297^ 








27 


30033 02 
5120 ^ 


2990637 o2 
102400 
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19899 /o 
3200 ^ 


257769 a 
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8739 


20 


10 


400 





9 
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45 
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9 
4 

















9 
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45 
8 














9 
4 



(5.17) 



in (|3.12p . which also rescales all operators, would lead to a matrix 



(5.18) 



The appearance of the /?-factors in ()5.17p leads to a "natural" form of Aij , according to the discus- 
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sion in section 13.31 as it ensures that the structure of the dominant solutions of appendix lA.ll stays 
intact, i.e., also after the redefinition the same combinations of [i and e"^ appear as before. 

Let us turn to the matrix Z from (j3.4p . The asymptotic solutions in the appendix have been 
chosen such that the submatrix of Z involving only the A < 4 operators vanishes identically. To 
obtain the other components, one would have to calculate more sub-leading terms in the dominant 
asymptotic solutions, but as stated in section [33} one can always choose a basis such that Zij = 0. 
This statement holds also after the operator redefinition given by (j5.17|) . Allowing for scheme 
dependence in the rotated basis, one would find from (|3.12p 



Calculating the two-point functions of the dual operators using ()3.8p is, of course, more involved 
and goes beyond the scope of this paper. It would be interesting to make the scaling with N^.^ ~ 
ln(A:/A) explicit, where A is the confinement scale, but certainly one would have to rely on numerical 
methods. We note that in |13[ [Tl] an approximation method was devised to determine the leading 
order term of the two-point functions in an expansion for large momenta. This method was also 
applied in [T6] . It would be interesting to see whether and how the renormalization procedure 
presented here would modify the results of |13l I14j . Another remark is that, in any case, the 
renormalization procedure devised here justifies a posteriori the pragmatic approach taken in [19] 
for the calculation of mass spectra in the KS system. Probably, it also applies to the calculations 
in |20l [2T] if it is possible to bring the system considered there into the form of a fake SUGRA 
system. 

It is also interesting to consider the counterterm matrix Uah- It is rather complicated, and we 
only give the leading terms in an expansion in e = e"^'^/'^. For Xij = 0, it is given by 




(5.19) 




(5.20) 



18 



with the submatrices 



UaxA — 



32 

" 15 
32 
5 



32 

" 5 
96 
5 



-g|o(32r2 + 148r- 873)/?2e2 



-g|5(32r2 + 148r- 873)/?2e2 
3/3e 



20 



20 i-^^ 



9/3e 
20 

20 '^'^ 

3(4r+17))3e 
16 



^3/2 



/ f(28r-31) if(2r + l) 



^73X3 = 



4x3 



1 



16 
3 



At + 5 



4r + 5 
/-|(4r + 17) 

8 
3 

\ 3 



Tf(2r+19) 



15' 

T|(2r + 19) 



W(r + 2) 



16 
9 



0{e)\ 
0{e) 



9/ 



(5.21) 



The entries of f73x4 and f74x3 lead to mixings between the fields in the glueball and gluinoball 
sectors. When considering non- vanishing Ajj, one notices that they are all scheme dependent. 

In general, the scheme dependent terms should only lead to finite contributions to the action. 
We have verified this explicitly for the A < 4 operators. Using the counterterm matrix with Ajj 7^ 
and considering non-vanishing sources only for the operators with A < 4, we find 



.(3) 



(5.22) 



with 



1 



^3l/3p4 200^ 

9 



( IOA45 — 4A44 2A45 ~ 2A54 + 5A55 

2A45 — 2A54 + 5A55 4A55 

IA47 + 5A65 — 2A64 — A46 IA57 — A56 + 2A65 

y 2A47 ~ 2A74 + 5A75 2A57 + 2A75 



IA47 + 5 A 



65 



2A 



64 



^46 



1^57 — A56 + 2A65 

5A67 — 2A66 
IA77 + 2A67 — A76 



2A47 ~ 2A74 + 5A75 
2A57 + 2A75 

A76 

4A77 



2 A77 + 2A67 



(5.23) 



These finite terms are analogous to the finite quartic counterterm in the GPPZ flow and the finite 
quadratic counterterm in the CB flow, cf. (j4.6|) and (j4.25p . respectively, after expanding them to 
quadratic order in the fluctuations. 

In addition to these finite terms there are also divergent contributions which are either linearly 
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diverging in e = e ^"^/^ or logarithmically. These are given by 
1. 



y(2) 



4-7 



_3l/3p4g2^0 X 

9 

/-|(r2 -3r + 5)/3 -|(4r-7)/3 



and 



(4r - 7)/3 








i(16r2 + 28T + 19) -4(2r + l) 

-16 




(5.24) 



with 



^(3,A=4) ^ ^2 



-L Z5_3 , 189^2 
16 """16 32 ' 



5355 I 38979 
128 256 



3 256t4+2368t3+4368t^-5664t-1821 
256 4t+1 



y(3,A=3) 



3 256t"^+2368t3+4368t^-5664t-1821 ' 
256 4t+1 

9 32r3+424r2+916T+371 
64 4t+1 



(5.25) 



/3^4 , 45,3 , 477,2 I 1959, 1239 4,3 , Ki _2 , 975, 83l\ 
y(3,A=3) = p + "2"^ + ~^ + ~8~^ ~ ~32~ + + ~^ ~ IF \ 

\ 4r3 + 51r2 + - ^ Ut"^ + 156r + 213 / ' 

Note that the linear divergencies are momentum independent for the A = 3 operators and pro- 
portional to for the A = 4 operators. Furthermore, the logarithmically divergent terms are 
proportional to k'^ for the A = 3 operators and proportional to for the A = 4 operators. All this 
is very reminiscent of the aAdS case, cf. sections 5.2 and 5.3 in [7]. There is, however, a difference 
in the fact that the logarithms appear in a much more complicated way, and they are even present 
in the linearly divergent terms. Although some of this may be an artifact of the choice of radial 
variable, this is consistent with the fact that the KS theory has no UV conformal fixed point. 

As mentioned above, all the entries of C/3X4 and U^xs are scheme dependent and thus only 
contribute finite terms to the renormalized action. This implies that one could have determined 
all the divergent terms for the glueball-sector and the gluinoball-sector separately. In other words, 
one can renormalize the KT theory without embedding it into the KS theory. This is plausible, 
as the KT background is a good approximation to the KS background in the asymptotic region, 
and the field theory divergencies are UV divergencies. This standpoint was also taken in |15l I16j . 
where the renormalization was performed just for the KT background. Indeed, it is easy to check 
that one can derive the diagonal components of the counterterms, i.e., U^xa and U^xs, by setting 
the last three components of the dominant solutions 01,03, 0.4 and 05 to zero when using (|3.3p . as 
well as the first four components of the dominant solutions 02, and dy. As explained at the end 
of section this corresponds to decoupling the glueball from the gluinoball sector. 

Finally, we would like to comment on the issue of VEVs. As we saw in the aAdS cases of 
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section 21 the response function of the background fluctuation W,j)/W encodes the VEV in those 
cases (up to an overall factor for the CB flow). We would like to see how this carries over to the 
case of KS. In order to derive the VEV from first principles, one would need the exact form of 
the counterterms linear in the fluctuations, which we have not determined yet. Thus, we can only 
take the cases of GPPZ and CB as encouraging examples and calculate, in analogy, the response 
coefficients of W^/W. It is straightforward to calculate 



W 



4r + 1 



-1 
1/3 
-2(r - 1/4) 

O3 



+ 



4e- 



4r + 1 



O4 

-4t + 1 
-4r+ l)(r- 1) 
\^ (t-2)(4t-1) J 



+ 0(e-2-) 



(5.26) 



Comparing this with the asymptotic solutions of the appendix we obtai: 

— = -2a5 - 407 + 2d6 . 



(5.27) 



This result suggests the interpretation that a combination of the two A = 3 operators has a VEV, 
which is in agreement with the field theory expectation of a condensate of the gluino bilinear 
[11^ I36j. However, this statement is again scheme dependent. The redefinition (|5.14|) leads to 



SiO-i + Xihi = Sja- + {xi - 5jXji)ai 



(5.28) 



and applying this to results in 



— = -2a5 + (-4 + 2A57)a7 + (2 + 2A56)a6 + 2A55O5 + 2A54a4 



(5.29) 



Let us apply the "naturalness" criterion on the form of the Ajj described in section 13.31 It would 
give A55, A54 ~ /3^, but /3 = in (j5.29p . so that the coefficients of 04 and 05 vanish. The coefficients 
of the dg and hj belonging to the A = 3 operators are more subtle, because the e""^ term in 05 
is independent of /?. On physical grounds we expect that there should be a natural scheme in 
which the VEVs for the A = 3 operators are not both vanishing simultaneously, cf. |1H I36j . It 
would be very interesting to understand how to determine such a preferred scheme, which might 
amount to extending the "naturalness" criterion of section 13.31 or to finding an equivalent of the 
supersymmetric scheme in the GPPZ flow. In this respect, it is interesting to notice that the 
term of 05 in (|A.5P can be written as e~'^/(4r + 1) x (04,-4,5,-9)-^ — 2d6. This suggests that 



^^In principle, there could also be a contribution from 03 in (|5.27|) . arising from the 0{e terms in (|5.26|l . which 
we did not specify in 05. However, this is scheme dependent, and one could always choose a scheme such that the 



coefficient of 03 in (|5.27|) vanishes. 
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the analog of the supersymmetric scheme in the GPPZ flow (which amounts to having a vanishing 
contribution of the sub-dominant solution to the dominant one, i.e. 02 = in (14.120 ) might be 
given by choosing A56 = 2 and A57 = in (j5.29p . We leave it for future work to make this argument 
more precise. Obviously, it would also be interesting to obtain the VEVs independently, using the 
linear terms of the action, but for this one would need the linear counterterms. 

6 Conclusions and Outlook 

In this paper, we made a first step towards the holographic renormalization of a general fake 
SUGRA theory, i.e., a theory whose scalar potential is of the form (|2.2p . One motivation comes 
from the aim to extend the program of holographic renormalization to backgrounds which are not 
asymptotically AdS, like the KS background. We devised a recipe how to calculate renormalized 
two-point functions for all the operators dual to the bulk scalars and discussed how the correlators 
depend on the renormalization scheme. 

Following the philosophy of [H |^ , our approach minimizes the effort in order to perform the 
renormalization of two-point functions. This means that we only determine the counterterms to 
quadratic order in fluctuations around a given background. Consequently, it is not obvious at the 
present stage how to obtain them from a fully covariant general expression upon expansion, if any 
such expression exists. Thus, our approach is very different from the one taken by [15\ [IB] in the 
context of the KT theory. Another difference is that we only consider the renormalization of field 
theories living on a flat space-time, i.e. we do not discuss any counterterms involving the space-time 
curvature. 

There is, obviously, much room for further development. First, it would be interesting to extend 
our analysis to other than two-point functions. We only made some preliminary observations 
concerning the VEV part of the one-point functions, based on the fact that the zero mode solution 
of the fluctuation equations seems to encode the VEV for the GPPZ and the CB flows. However, 
a more systematic understanding involving an action and counterterms linear in the fluctuations 
would be worthwhile. A full understanding of the VEVs would also require a way how to choose a 
"natural" scheme, for example by extending the naturalness criterion of section 13.31 or by finding 
other selection criteria, similar to demanding supersymmetry in the case of GPPZ. Moreover, it 
would be interesting to describe also the non-linear interactions of the gauge-invariant variables 
given in [26], including the metric fluctuations, in the effective action approach used in this paper. 
A crucial check of consistency, which we have not touched upon, is the verification of Ward identities. 
At first sight, it appears that they are inherent in the use of the gauge-invariant variables, but a 
more detailed analysis might also shed light on the counterterms of other than second order in the 
fluctuations. 

Second, it would be desirable to scrutinize the assumptions made and look for settings where 
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they do not hold. A particular candidate would be the Maldacena-Nunez system [37], for which 
the consistent truncation to a fake SUGRA system has been considered in [25\ |26] . It would also 
be interesting to see whether and, if yes, how the systems considered in [10] could be described in 
our approach. 

Third, it would be interesting to consider non-Poincare sliced backgrounds and to include fields 
other than scalars, e.g., vectors. We hope to come back to some of these questions in the future. 
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A Asymptotic KS Solutions 

A.l Dominant Solutions 

The dominant asymptotic solutions, up to order e"^"^/^ relative to the leading term, are (in mo- 
mentum space) 
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